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ABSTRACT 

We derive the distribution of interferometric visibility for a source exhibiting strong diffractive scintil- 
^^ . lation, with particular attention to spectral resolution at or near the Nyquist limit. We also account 

7-H ' for arbitrary temporal averaging, intrinsic variability within the averaging time, and the possibility of 

f^ \ spatially-extended source emission. We demonstrate that the interplay between scintillation and self- 

^SJ . noise induces several remarkable features, such as a broad "skirt" in the visibility distribution. Our 

» , ' results facilitate the interpretation of interferometric observations of pulsars at meter and decimeter 

^ il wavelengths. 

■^ ■ Keywords: methods: data analysis - methods: statistical - pulsars: general - scattering - techniques: 

high angular resolution - techniques: interferometric 

(N ; 

^ ; 1. INTRODUCTION 

HH \ Radio astronomy is now exploring two limiting regimes. First, with the construction of ever-larger collecting areas, 

f-H . noise intrinsic to the source (i.e. self-noise) can be a major, or even dominant, component of the noise budget. Second, 

f-^ > observing systems at decimeter wavelengths can now record baseband data over wide bandwidths and with sufficient 

I I bit-rates to effectively capture the full information content of the signal (up to the fundamental limit described by 

O , the Nyquist sampling theorem). These advances motivate mathematical descriptions of the new regimes and the 

^ • development of novel techniques that optimally utilize the signal inform atio n. 

c/3 ' These considerations prompted a previous work (Johnson & Gwinn 2012. hereafter JG12), in which we derived the 

I ^1 , probability density function (PDF) of flux density, sampled at the Nyquist rate, for a scintillating source. We accounted 

for background noisqj, arbitrary temporal averaging, the possibility of decorrelation of the scintillation pattern within 

T— I I the averaging time, and spatially-extended source emission. We also outlined tests to identify self-noise, as well as 

^ . a mechanism to detect rapid intrinsic variability of a signal. We then demonstrated the success of this description 

CO on 800 MHz observations o f the Vela pulsar and, thereby, obtained a spatial resolution of approximately 4 km at the 

<^i pulsar (jJohnson et al.|[2012[ ). 

OO , We now extend these mathematical results and techniques to interferometric visibility, which preserves much of the 

^^ ' statistical framework. This extension involves a pair of stations, which observe independent background noise but 

^^ ■ correlated scattering kernels and identical source noise. However, these mild modifications lead to substantial changes 

^^ in the analytical results for the distribution functions. Furthermore, the interferometric visibility is complex, and thus 

^f^ I is characterized by a two-dimensional PDF, whereas the flux density is a positive real number. We frequently utilize 

^-H ' projections as a natural mechanism for exploring the transition between these domains. We also offer several other 

l" , close analogs with common metrics for flux-density, and we provide tools to quantify the complex noise in sampled 

. !^ ' visibilities. 

K^ Pulsars are some of the richest targets for interferometry, and we have tailored our treatment to account for their 

5h ' extraordinary variability. In particular, we emphasize results for single-pulse studies with Nyquist-limited resolution. 

C^ . We then connect to traditional limits through asymptotic forms of our results, as the averaging is increased and self- 
noise becomes negligible. However, our principal goal is to facilitate precision tests using interferometric observations, 
particularly of pulsars, that can robustly distinguish between intrinsic and extrinsic characteristics and can sensitively 
probe delicate modifications of either. 

1.1. Assumptions and Strategy for Comparison with Observations 

Our scope is deliberately broad, intended to encompass the majority of pulsar observations at meter and decimeter 
wavelengths. However, to assess the suitability for comparison with any particular observation, we now outline the 
assumptions underlying our results. Our physical requirements are that 

• The source emits amplitude-modulated noise ()Rickettlll975[ ). 

• The scattering is strong. That is, the diffractive scale r^ = X/9 is much smaller than the Fresnel scale rp = vAD, 
where A is the observing waveleng th, 9 is the angular siz e of the scattering disk, and D is the characteristic 
distance to the scattering material (jCohen fc Cronvnlll974[ ). 

I michaeltdh@physics . ucsb . edu, cgwinn@physics . ucsb . edu l 

^ By background noise, we mean the sum of all noise other than that of the source (e.g., receiver, spillover, atmosphere, and background 
sky noise). 
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The first assumption is quite robust for most, if not all, astrophysical sources, simply because the superposition of 
many independent radiators will produce a signal of this form. Also, the majority of meter and decimeter pulsar 
observations fall easily within the regime of strong scattering: r^/r-p <^ 1. 
On the other hand, the instrumental assumptions that we require are that 

• The data explore a large representation of the full ensemble of diffractive scintillation: {B / /S.v^j'x (iobs/Atd) S> 1, 
where {_B, iobs} are the total observational bandwidth and duration, and {Ai'd, Atd} are the characteristic 
bandwidth and timescale of the scintillation pattern. If B / /S.va ^ 1, then intrinsic single-pulse amplitudes may 
be estimated, allowing a comparison between the data and our models that requires no fitted parameters. 

• The data are coherently dedispersed, baseband shifted, and "fringed" with an appropriate phase model. 

• The time to form each spectrum, or the accumulation time iaco is much longer than the pulse-broadening 
timescale ip but much shorter than the characteristic scintillation timescale Atj- 

The last assumption is perhaps the only one that is atypical of modern observations and processing. However, it 
ensures that the spect ra reflect the convolution action of the scattering on the in trinsic pulsar signal with a stochastic 
"propagation kernel" (jHankinsI 119711 : IWilliamsonI 119721 : IGwinn fc Johnso"nll2011[ ). For pulsars, a particularly elegant 
limiting case is the formation of spectra that include all the pulsed power. This convolution leads to extremely 
general results; we do not assume thin-screen scattering or a particular spectrum for the density inhomogeneities of 
the interstellar medium (ISM), for example. 

Our work also addresses the physical possibility of an extended emission region and the instrumental possibility of 
averaging (or integration) of the calculated visibilities. For the former, we assume that the source is a small fraction of 
the magnified diffractive scale (see Sj3|); for the latter, we assume that the averaging is of N ind ependent cross-spectra 
say, f rom different pulses), over a timescale much shorter than Aid (i-e. the "snapshot image" of lNaravan fc GoodmanI 
1989( )). While frequency averaging is analogous, it inherits additional information from non-stationary signals and 
requires the specification of individual pulse profil es, wh ereas a description of temporal averaging requires only the 
phase-averaged source intensity for each pulse (see IJG12I for details). Hereafter, we exclusively use N to denote this 
degree of temporal averaging. 

1.2. Relation to Previous Work 



For non-scinti ll ating sources, the statistics of interferometric visibility are well-known; see Eoran (197 61) or 
[Thompson et al.l (|2001| ). In particular, the self- noise has been carefully charac terized. For example, Kulkarni. (jl989f ) 



analyzed the noise in synthesis imaging for sources of arbitrary strength, while lAnantharamaiah et al.l (|1989| ) studied 
the noise, with an emphasis on extremely bright sources. 

Scintillation complicates th e sta tistics, and even the no i se-free (i.e. infinite-averaging) case warrants careful attention. 
iNaravan fc GoodmanI ([198!/) and 'Good man fc NaravanI (11989*) analyzed this li mit via both numerical and analytical 
techniques, with an emphasis on moments of the distribution of visibility. Also, iGwinnI (|2001[ ) calculated the full PDF 
of interferometric visibility for a scintillating source in this limit and accounted for the effects of an extended emission 
r egion. 

IGwinn et al.l ()2012af) and IGwinn et al.l ()2012b[ ) then incorporated the contribution of both self- noise and background 
noise for modest averaging {N > 20) by applying the central limit theorem to the averaged visibilities. The present 
work extends these ideas by adopting a different strategy, following JG12- Namely, we establish the statistics for 
individual spectral samples and then determine the effects of averaging by a convolution of independently drawn 
samples. This method permits the treatment of data with arbitrary averaging, iV, including no averaging (i.e. N = 1), 
and requires no assumptions about the nature of the intrinsic variability or of the scattering material. 

To put our results in the context of these earlier efforts, we also present asymptotic forms of our equations for strong 
signals and high degrees of averaging. However, our principal goal is to facilitate direct comparisons with observations. 

1.3. Outline of Paper 

In ^ we focus on visibility statistics for a scintillating point source. We account for temporal averaging, but assume 
that the scintillation pattern is fixed within each average. Under these constraints, we first outline the essential 
statistical framework for the visibility statistics (i 32.ip and calculate moments and noise of the visibility statistics within 
a fixed scintillation element (i i2.2|) . We then outline approximation schemes for the PDF of visibility during a fixed 
scintillation element ( §2.3|) and demonstrate how to calculate the PDF of visibility after including the scintillation 
ensemble ( ^2.4p . We give examples of the visibility statistics in limiting regimes ( ii2.5p and demonstrate that the 
combination of self-noise and scintillation introduces a "skirt" in the visibility PDF that dominates its asymptotic 
form, regardless of the source strength or the baseline ( ^2.61) . We also provide a prescription for estimating the 
self-noise (^Ji). 

Next, in SjSj we quantify the infiuence of a spatially-extended emission region on this PDF and discuss the potential 
of interferometry to resolve the size and anisotropy of such emission. We demonstrate that the effect of a small 
emission region depends only on two parameters, regardless of the scattering geometry; however, the translation to a 
dimensionful size at the source depends on the scattering geometry. We also derive a simplified version of the PDF of 
visibility for the special case of a zero-baseline interferometer. 

Finally, in 21 '^^ summarize our results and outline some observational prospects. 
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2. PDF OF VISIBILITY 



We now derive the expected PDF of interferometric visibility arising from a scintillating point source. We allow 
arbitrary temporal averaging, in sets of N independent cross-spectra, but assume that the scintillation pattern is fixed 
during each average. All of our results apply to scalar electric fields (i.e. a single linear or circular polarization) that 
conform to the assump tions d escribed in m.l\ 

Our notation follows IJG121 Namely, Zx denotes a circular complex Gaussian random variable with unit variance, 
indexed by x; Gx denotes an exponential random variable with unit scale. We use Pi) to generically denote a PDF 
with respect to the given variables and parameters, and we present PDFs P{'w) of complex quantities w with respect 
to the metric dRe[u']dIm[u;]. Occasionally, we employ the shorthand w^ = Re(u'), Wi = lm{w). 

Furthermore, to help visualize the visibility PDF, we rely on projections. These projections are most enlightening 
along the real axis K- = R-e(F), where they emphasize the relative effects of differing baselines. We favor the following 
two projections: 



P{y,-N) = j dV,P{V-N), (1) 

Q{V,-N)= [dViV^P{V;N). 



The first projection identifies the concentration of density toward greater real part, which reflects the average visibility; 
the second (weighted) projection quantifies the imaginary spread of density, which reflects the influence of scintillation 
and the noise from both the background and the source. 

2.1. Field Statistics 

In §3.1 of lJG12l . we derived the electric-field statistics for a scintillating point source. These statistics depend on the 
amplitude-modulated noise of the source, the strong scattering of the ISM, and the receiver noise. We now review the 
basic ingredients of this description and establish the necessary notation. 

For the electric field, the amplitude- modulated noise takes the form y/AjIgfiei. Here, Aj is a dimensionless amplitude 
factor, indexed by pulse, that accounts for pulse-to-pulse variations. Is is a constant characteristic scale of source 
intensity, fi is a power-preserving envelope, and Ci is white Gaussian noise of unit variance. Thus, the pulse profile is 
simply l/ip, the gated signal has mean amplitude A jig, and Ci accounts for the noiselike nature of the emission. Note 
that this treatment accommodates arbitrary variability of the pulsar, such a s the possibility of correlated pulse-to- 
pulse variations, log-normal amplitude statistics, or nanosecond-scale bursts (jRickett et al.l[l975t iCairns et al.l 120011 : 
iKramer et al. 2002 ). 

During a fixed scintillation element, the scattering acts to convolve this intrinsic emission with a stochastic prop- 
agation kernel: girji. In the strong-scattering limit, the form of this kernel is similar to the emission of the pulsar. 
Namely, a power-preserving envelope gi modulates Gaussian noise r]i. This envelope is more commonly described by 
its squared norm: the pulse-broadening function \gi\'^. 

Finally, an observer samples the propagated signal in the presence of white background noise: ^/TZPi- Here, /n 
is a constant, characteristic scale of the background noise, and (3i is white Gaussian noise of unit variance. If the 
background noise changes significantly with time, then a changing scale may be added, analogous to Aj. 

The observed scalar electric-field time series Xi and its Fourier-conjugate spectrum Xi are thus given by 

x^ = ^/A~h [(/e) * {gTl\ + ^/I^|3^ ^ i^ ^ ^/A~Ts (/*?)_ (g * ??), + V7n/3z, (2) 

where a tilde denotes a Fourier conjugate variable. 

Because q, ry^, and /3i are mutually independent (circular complex Gaussian) white noise, a single spectral sample 
Xi is of the form yM^TgZfZg -|- ^/T^Zi^, where Z[ = (/ * e)i, Zg = {g* fj)i, and Zb = A are each circular complex Gaussian 
random variables with unit variance. If the scintillation is held fixed (i.e. Zg = const.), then the intensity \xi\'^ is drawn 
from an exponential distribution with scale Ij = Ajls\zg\'^ + /„ . 

For interferometric visibility, the observer measures the cov ariance of the e lectric fields at two stations, which we 
denote by unprimed and primed variables: Vi = Xix'* (Thomp son et al.ll200lD . The electric fields at the two stations 

arise from identical intrinsic emission, f * e' — f * e; however, a difference in sensitivity or gain between the stations 
will affect the overall scale: Ig ^ I^. The background noise is independent ((zbz(j*) — 0) and with different variance 
{In ^ /„) at the two stations. Finally, the stochastic part of the propagation will also differ, with a baseline-dependent 
correlation pg = l^Zgz'*). 
Hence, the average of N visibilities from different pulses takes the form 

1 ^ 
V^ = ^ XI (V'^Vs^f J"^g + v/^^bj- j [y/AjI'^zf^jz'^ + ^/T'nZ^yj . (3) 

i=i 

Even in the zero-baseline limit (z' ~ Zg), Eq.[3]differs from the corresponding intensity result because of the assumption 
of independent background noise at the stations. 
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As Eq. 13] shows, the measured visibihty is the average of N random variables, each of which is the product of 



two correlated circular complex Gaussian random variables with respective variances Ij 
Ajl!^\zl\'^ + /,j and correlation 



Ajlslzgl' 



and Ij 



f g J 




(4) 



This correlation changes with scintillation and is complex because of covariance between the real and imaginary part 
of electric fields at the two stations. 

For pulsar observations, one can examine the mean of the off-pulse and on-pulse spectra at each station to estimate 
the parameters Jg, /g, /„, /^, and Aj. These measurements fully characterize the visibility statistics of Eq. [31 

2.2. Moments and Noise of Snapshot Visibilities 

Using Eq. [3J we can evaluate moments of the visibility distribution for a "snapshot image" (i.e. the scintillation 
variables Zg and z' are held fixed). These moments include the effects of pulsar variability and self- noise. For example. 



J T'z z'* 



(Re(y)2) 
(Imivr) 



{V) = {A)m 



(5) 
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These expressions present ensemble averages over the noise of the pulsar and background, while the set of N pulse 
amplitudes {Aj} and the scintillation factors are held fixed; (. . .) denotes the average over noise, whereas {. . .)n denotes 
the average over the N pulse amplitudes. 
We can also quantify the noise using these moments. Within a single scintillation element, the noise takes the form 



im')^ 




(6) 



/f/n)|W| + l|l 



{6A^) 



N 



{A)% 



\{V)f 



= 2b„ + 2bi\{V)\+b2\{V)f 



Thenoise is a quadratic function of the signal; hence, pulsar field statistics are heteroscedastic ()Gwinn et al.ll20lTl 
l2012al). In part icular, this equation demonstrates the contribution of the source to the noise, as originally investigated 
bv lDickj (|1946D and described by the familiar radiometer equation. Observe that, for the special case of a zero-baseline 
interferometer, the noise coefficients, &i, are independent of the particular scintillation element {zg^z'}. 
We can similarly obtain the variances parallel and perpendicular to the mean visibility: 



SV^f)^bo + bi\{V)\+b2\{V)\\ 
{5Vl)=bo + bi\{V)\. 



(7) 



Thus, the noise scales quadratically in phase with the signal but only scales linearly at quadrature to the signal. Figure 
[T] illustrates this noise behavior. 

Because they neatly separate the contributions of the background and source, the noise coefficients provide a valuable 
mechanism to study and quantify the various types of noise. The self- noise coefficient, 62, is particularly useful because 
it contai ns information about the intrinsic variability, which does not affect the aver age spectrum or correlation 
function (IGwinn fc Johnsonll201lD . Other types of noise, such as quantization noise (|Colelll968: [Thompson et al.l 
1200 It IJenet fc Andersonlll998l ). will modify the coefficients but preserve the quadratic form (fGwinnI 120061 ). In i }^771 we 
describe a procedure to estimate the self-noise. 

2.3. Approximating the Visibility Statistics within a Scintillation Snapshot 

We now derive the PDF of visibility for samples that are collected and averaged within a single scintillation snapshot, 
characterized by {zg,z'}. Because the field statistics take a rather general form, P{V\zg,z') is simply the iV-fold 
convolution of the distribution of the product of correlated complex Gaussian random variables; we relegate the 
derivation and details of this distribution to the appendix (in particular, ! jA.l[) . 

We also present several approximation strategies, because the convolution of visibilities that are not statistically 
identical has no convenient analytical form. These strategies provide accuracy that is sufficient for most applications 
and constitute the analytical foundation for our subsequent results. Moreover, the first approximation that we derive, 
the i.i.d. approximation, is exact for A^ = 1. 
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Figure 1. Characteristics of complex noise SV for snapshot visibihties. (ie/i) EUipses denote the standard deviation of the Gaussian 
distribution of noise, centered on the mean visibility. The major axis is or iented with the phase of the visibility and scales quadratically with 
the signal, while the minor axis scales linearly with the signal (see il2.2|l . The plotted noise corresponds to noise on a short baseline with 
all source and background intensities unity, and N = 100. (right) Average complex noise {5V'^{V)), estimated as described in ^2.71 The 
increasing vector length with \V\ shows the increasing influence of self-noise, while the change in direction reflects the varying orientation 
of the noise ellipse. 

2.3.1. The i.i.d. Approximation 

We derive our first approximation by assuming that the averaged visibihties for each set of N pulses are independent 
and identically distributed (i.i.d.). To achieve this condition, we treat the pulse amplitudes as constant for each set 
of N averaged visibilities: Aj i— > ^ = {A)i^. Nevertheless, this approximation preserves some information about the 
pulsar variability; the i.i.d. approximation is exact when iV = 1, for instance. 

After this replacement, Ea. lA3l gives the PDF of visibility: 



,N 



p{y-N\z^,z') 



p = z^z^A\jjj^, 



2^'.(iV-l)! a^+i '^' ^^-' r — j ""P r ^T- 



'IT, J=A/,|zgp + /„, T = AII\z:Y + I[^ 



(8) 



In this expre ssion, K^ {x) is the modified Besscl function of the second kind. This approximation and its equivalent 
for intensity pG12l §3.2.1) provide convenient tools for analytic work. 

2.3.2. The Gaussian Approximation 

As the number of averaged samples N -^ oo within a fixed scintillation pattern, the PDF of visibility approaches 
an elliptical complex Gaussian distribution with mean V ~ (A) n yH^Jlzgz'* and variances determined by Eq. [71 
iGwinn et al.l (|2012aD and lGwinn et al.l (|2012bD described this limit and result, and verified the noise prescription using 
observations of the Vela pulsar. 

In general, the linear term, bi, of the noise polynomial depends weakly on the scintillation {zg,z'}; however, for a 
short baseline, the noise ellipse depends only on the mean visibility. 

2.4. The PDF of Visibility 

We now derive the PDF of visibility, when the data explore a representative ensemble of the diffractive scintillation. 
This result relies on both the PDF of visibility within each scintillation snapshot (i i2.3p and the PDF of the scintillation 
random variables {zg, z' }: 



P{V- N)= j P(V- iV|zg, z;)F(zg, z'g)(fzg<fz'g. 



(9) 



Now, Zg and z' are drawn from a distribution of circular complex Gaussian random variables with some correlation 
Pg = (zgz'*)s. Note that we use the subscripted brackets (. . .)s to designate an ensemble average over the scintillation. 
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Figure 2. Visibility projections as a function of baseline and signal-to-noise (SNR) for Af = 1. We have set the pulse amplitude A to be 
unity, the receivers to have identical signal-to-noise ratios (/s = /^ and /n = ^4), and the single-dish intensities to be fixed: 7s -|- 7n = 1- 
The source contribution can increase the variance about the real axis on long baselines as a result of scintillation. This combination of 
projections demonstrates that, even without averaging, a strong scintillating source may be readily detected regardless of baseline. The 
short, medium, and long baselines correspond to pg = 1.0, 0.5, and 0.0, respectively. 



Although a single realization of the scintillation pattern has a complex mean visibility z^z'* we assume that the 
ensemble-averaged mean visibility pg = ^Zgz'*) is real. In practice, this assumption merely reflects an appropriate 
added, constant phase to the delay model. This mean visibility depends on the baseline, scattered image size, and 
observing wavelength, and can be expressed in ter ms of the phase structure function D^{\i) of the scattering medium 
(|Tatarskiilll97ll : ILee fc JokipiJ|[l97l IRickettlUOOOl) : 



Pg = exp 



-2^0(b) 



(10) 



The joint PDF P{zg, z' ), the standardized bivariate complex Gaussian distribution, follows easily from a distribution 
of four correlated (real) Gaussian random variables that correspond to the real and imaginary parts of Zg and z' ■ see 
iGoodmaiil (|1963[ ) or iGwiniil ()200lD . for example. In terms of the scintillation norms, r = |zg| and r' = \z' \, and their 
relative phase 9 = arg (zgz'*) , we have 

P{r,r',0) = —— :7V "^xp 



TT (l 



Pi: 



(r^ + r'2 - 2pg rr' cos 9) 



(11) 



2.4.1. The i.i.d. Approximation 

For the i.i.d. approximation f ^2.3.1l) . we can further reduce the expression for P(y;N) by integrating 9 because 
arg(p) is the only quantity in Eq. [S]that depends on 9. We then obtain 
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In this expression, /q (x) is the modified Bessel function of the first kind. The visibility projections and their general- 
izations, such as the projected imaginary variance J dViV{^P{V), similarly follow from the results of ! jA.2l Figure [2] 
illustrates the effects of baseline and signal-to-noise on the projections. 
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Figure 3. Visibility projections P{Vr) and Q(V|.) for N = 2, 10, and 50, and the corresponding errors 5P and 5Q in the Gaussian 
approximation. We have set the signal-to-noise ratio to unity (Is = /g = /n = /^ = 1) and assigned a moderate baseline (pg = 0.5). We 
have set all pulse amplitudes to be equal to one, so the i.i.d. approximation is exact. 



2.4.2. The Gaussian Approximation 

We now include the effects of scintillation for the Gaussian approximation (i j2.3.2l) . The distribution of v isibility, 
after including the scintillation ensemble but before including the effects of noise, can be written (see Eq. lAip 



P{V;N ~>oo) 



1 



-Ko 
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exp 



2Pg 



MV]\ 



(13) 



where n = {A)N^yl^ and pg = (zgz'*). IGwinnI (|2001l ) described this hmit and result. 

As noted in ^2.3.21 the noise only depends on the mean visibility when the baseline is short (relative to the difFractive 
scale). Hence, the Gaussian approximation, including the scintillation ensemble, becomes 



P{V; N) « / P{V'; N ^ oo)Pnoisc("t^ - V , V')d?V' . 



(14) 

In this case, Pnoisc(V^, Vb) denotes the elliptical Gaussian distribution of noise centered on Vq. The orientation of the 
ellipse is given by the phase of Vq; the major and minor axes depend on both Vq and the noise coefficients {&0j ^i; ^2}, 
as derived in ij2.2l 

From Eq.[T31 we see that the visibility projections after including the effects of noise can be applied directly to Pnoiso- 
For example, the real projection is given by 



P{V,]N)fv / P{V']N ^00) 



1 



1 



/2i^ho + bi\V'\+h2 



: exp 



(Vr~Vjf 



2{bo + bi\V'\+b2\V'\'^cos^, 



d?V\ (15) 



where </> = arg V . iGwinn et al.l ()2012aD and lGwinn et al.l (|2012b[ ) used these representations to characterize the PDF 
of visibility for observations of the Vela pulsar. 

The Gaussian approximation is also effective on longer baselines, although the noise depends on the scintillation pair 
{zg, z' } rather than just the mean visibility. However, in this case, we simply replace 61 by its value when Zg = z' ■ in 
this case, the noise only depends on the mean visibility, and Eg. 1141 applies. Figure |3] shows the error in the Gaussian 
approximation for various degrees of averaging on a moderate baseline. 

2.5. Examples 

Because the form of the visibility PDF is rather opaque, we now refine it for several cases of interest. The zero- 
baseline interferometer (t ^2.5.ip offers the most substantial reduction of complexity and relates neatly to the analogous 
results for intensity. At the opposite extreme, an infinite baseline ( ii2.5.2[) also leads to a greatly simplified form, albeit 
with no fundamental decrease in numerical difficulty. Finally, the regime of high signal-to-noise ( ^2.5.3p highlights the 
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influence of self- noise on the PDF and lays the foundation for calculating the asymptotic (|T^| — >■ cx)) dependence (see 
i j2.6p . For simplicity, we derive all results of this section using the i.i.d. approximation. 



2.5.1. Zero-Baseline 

The zero-baseline limit for visibility is particularly simple: pg — >■ 1. In this case, Zg = z' and the effects of scattering 
depend only on the scintillation "gain" G = \zg\'^: 
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where P{G) ~ e^ . This representation has two distinct advantages: it requires a single integral, and it can be easily 
modified to any alternate distribution of scintillation gain, as might occur from an extended emission region or weak 
scattering, for example. 

We can also use the results of ijA.2l to evaluate projections of the visibility distribution without requiring an additional 
numerical integral: 
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In these expressions, k^ix) is the modified spherical Bessel function of the second kind (jArfken &: Webeill2005[) . 
As the averaging increases, the distribution of visibility approaches P{G): 

P{V- N^oo) = -^e-^°[^l/<^>0(Re[y])(5(Im[l^]), 
where 9{x) is the Heaviside function, and S{x) is the Dirac delta function. 



(18) 



2.5.2. Long-Baseline 

On a sufficiently long baseline, the interferometer completely resolves the scattering disk, and the respective prop- 
agation kernels at the two stations become completely independent: Pg — >■ 0. Although ground-based VLBI can only 
marginally achieve this regime for a few of the most heavily scattered pulsars at meter and decimeter wavel engths (e.g. 
Gwinn et al.lll993[) . space VLBI with RadioAstron can easily resolve the scattering disks of many pulsars (|Kardashevl 



2009|) . A precise understanding of the visibility statistics for ultra-long baselines will help to maximize the information 



that can be gleaned from space VLBI. 
In this limit. Eg. [T2] becomes 
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Observe that this distribution depends only on \V\, as necessitated by the phase invariance of this limit. Although the 
mean visibility is zero, the effects of scintillation can still be substantial, as Figure [2] illustrates. 

2.5.3. High Signal-to-Noise 

In the limit of infinite signal-to-noise, the fields at two stations will exhibit identical noise, multiplied by their respec- 
tive complex scintillation gains. In terms of the field statistics (i i2.ip . each visibility is the average of N exponential 
random variables with means Aj^IsI^Zgz'* . 

For simplicity, we apply the i.i.d. approximation. In this case, each averaged visibility is drawn from an Erlang 
distribution, multiplied by the scintillation gain, which rotates the phase to arg (zgz'*). Expressed in terms of the 

magnitude and phase of the visibility, V = \V\e'^'^, this distribution is 
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The definitions of r, r', and 9 are equivalent to those in ij2.3.1l When combined with the scintillation ensemble, 9 can 
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be trivially integrated against its delta function to give 
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To obtain the second line, we transformed to polar coordinates {r = icosip^ r' = Isinilj} and integrated over £. Thus, 
as for the zero-baseline case, we have reduced the integration to a single dimension. We next use this form to determine 
the asymptotic behavior of P{V). 

2.6. Asymptotic Behavior 

We now derive the behavior of P{V; N) as |F| — )• oo. To proceed, we first consider the high signal-to-noise results. 
At large \V\, the Bessel function of Eg. [2T] approaches an e xponential, so we can apply the method of steepest descent 
to approximate the integral over ijj (jArfken fc Weben 120051 ): 
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Here, the constant of proportionality is also a function of (j). 

On the other hand, the distribution of visibility for purely background noise has the asymptotic form (see Eq. IA3 
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Thus, the stronger, scintillation-induced "skirt" of source power will dominate the PDF of visibility in the asymptotic 
regime defined by 
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Of course, the projections also refiect this remarkable behavior. For example, the projection onto V^ = Re(T^) follows 
by applying the method of steepest descent to these asymptotic forms. We thereby obtain 
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Once again, a scintillating source extends the wings and introduces a baseline-dependent asymmetry. In this case, the 
constant of proportionality is a function of sign(T/r), and the asymptotic regime is determined by 
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(26) 



Figure 0] demonstrates this asymptotic behavior of the projections on various baselines. 

Observe that both these asymptotic forms have poles for pg = ±1. This behavior reflects the fact that, in these 
cases, the signal power is restricted to the real half-line that matches the sign of the mean visibility, so the asymptotic 
behavior elsewhere corresponds to that of pure background noise. 

Perhaps most surprisingly, in the limit N ^>- oo (given by Ea.ll3p. the asymptotic behavior is exponential, regardless 
of the signal-to-noise. Hence, the broad exp(— -^/JV^) skirt arises from the delicate interplay of the scintillation and 
the self-noise. 

2.7. Estimation of the Self- Noise 

As a last application of these visibility statistics, we now present a prescription that can be applied to quantify 
self-noise. As we have discussed, self-noise provides a powerful diagnostic of a signal, particularly when characterizing 
intrinsic variability. For example, source variability on timescales shorter than tacc induces correlations in the spectral 



10 



Johnson & Gwinn 



1 
10-1 
10-2 

10-3 

-i 10-4 
> 

oT 10-5 

10-6 
10-^ 
10-B 



eI 1 1 1 1 1 

r .^ 


|iiii 


" -'// 


'-XV : 


— '' / 


\x>\^ ^ 






- .-'' / / 


-^^ 


- --'' ^ / 


\ z 


- .-' / / 


"-.-! 


-' ^ / 




5- / / 


-5 


: / / 


Baseline \ 


/ / 






Short I 


=- / 


Medium ~| 


r / 


Long n 


"i 1 F 1 1 1 1 


1 1 1 1 1 1 1 r 




10 




V 



10 



I I I I I /I I I I I I I I I I I I I I I I 



-20 



10 




V. 



10 



20 



Figure 4. Exact projections P(V'r) and Q(Vr) for A^ = 1; the source and background amplitudes are unity. The broad scintillation-induced 
wings are distinctive at large IVrl- However, for zero-baseline and negative Vr, there is only the ex pone n tiall y- falling noise contribution, 
reflecting the absence of source power on that half-line and, mathematically, the singularities in Eqs. 1221 &: 1251 



noise, without modifying the mean spectrum (jGwinn fc Johnsodl20lil ). We now outhne one technique that estimates 
self-noise and thereby facilitates a detection of such corre lations. 

Our method is similar to its analog for intensity pG12l . §5.2). Namely, we use pairs of nearby samples to calculate 
finite estimates of the signal and noise. Each pair of samples is assumed to be within a single scintillation element. 
We will assume that the pair consists of samples with uncorrelated self-noise (e.g. pairs from different pulses, or from 
the same pulse with negligible intrinsic modulation). We calculate the noise in pairs of samples as a function of their 
mean. However, because the visibility is complex, the noise must be treated as a vector quantity as it has different 
behavior in phase and at quadrature with the signal (see §2.21 and Figure [1]) . We find it advantageous to work with 
analytical expressions of the involved complex quantities and, thus, define 



bV'^{V) = 2 
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We again work within the i.i.d. approximation. The distribution of visibility is then given by Eq. 

5V^{V) 



(27) 



31 and we obtain. 



iV- 



(28) 



As for intensity, this method for estimating the self-noise agrees with the exact expression with iV — > iV + 1/2. 

If the averaged visibilities are from different pulses, then pulse-to-pulse variations can contribute additional noise. 
If the averaged visibilities are from the same pulses, then intrinsic variations on timescales shorter than iacc induce 
correlations in self-noise, and th ereby decrease the m easured noise; if A^ > 1 then pulse-to-pulse variations within the 
averaging will increase the noise. I Gwinn et al.l (|2011[ ) applied similar tests to infer short-timescale variability (< 300^s) 
of PSR B0834-I-06. 

3. EFFECTS OF AN EXTENDED EMISSION REGION ON THE VISIBILITY PDF 

We now derive the modifications to the PDF of visibility from a spatially-extended emission region. Such emission su- 
perimposes many slightly-offset copies of the diffraction pattern at the observer, suppressing the observed scintillation. 
Optical scintillation provides a familiar demonstration: "Stars twinkle, but planets do not." Emission extending over 
a region much larger than the diffractive scale r^ (see ijl.ll) quenches the scintillation, and so difi'ractive scintillation 
studies can effectively probe emission scales < r^. 

3.1. The Effects of an Extended Emission Region on the Field Statistics 

If the emission spans a transverse size ^rd, then the modification to the field statistics can be derived quite generally. 
Explicitly, in terms of transverse source coordinates s and the notation of W2,.l\ the observed electric field takes the 
form 

(29) 
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We set the origin of the coordinates s so that J cPs s^(s) = 0. Because the source intensity /(s) is assumed to be 
confined within a region <^ rj, Zg(s) wih only vary shghtly and we may expand to hnear order: Zg{s) sa Zg(0) + 

The source term in Eq. [291is then a convolution of three complex Gaussian random variables. Because of our choice 
of origin for s, at linear order these three random variables are mutually uncorrelated during a fixed scintillation 
pattern. In addition, the scintillation random variable, Zg(s), is uncorrelated with its spatial derivatives, so the scales 
of the three respective variances are also mutually independent at linear order. Combining these characteristics for 
the pair of receivers then gives the form of the visibility field statistics to leading order: 



iV 1 + 7s4 + 7s_2 ^ '- 



(30) 



^/A~T^(zfjZg + ^TlTZf^ijZgi + y/%^Zf^2jz'g,2) + V^4„ 



Here, the paired scintillation random variables are correlated, pg 



and 



\^g,i^g,ii 



, but all other pairs of 



random variables {zx,Zy} are uncorrelated. We have chosen the scaling prefactor so that the intensity measured by 
either receiver is unaffected by the extent of the emission. Our expansion parameters, the dimensionless subsidiary 
scales 7s, i ^ 1, contain information about the transverse extent of the source emission. More specifically, these scales 
are proportional to the spatial standard deviation of integrated flux density. For example, spatially-offset, pointlike 
emission sites need only to emit within the same accumulation time, but not necessarily at the same retarded time, to 
affect the scintillation statistics. Hence, even for emission that is highly beamed a nd, thus, instan t aneou sly pointlike, 
this method can identify a transverse size that relates to the emission altitude; see iJohnson et all ()2012D . 

Within a fixed scintillation pattern, Eq. [30] is the A^-fold convolution of products of complex circular Gaussian 
random variables, as for a point source. Thus, applying the i.i.d. approximation, we see that the PDF of visibility 
within each scintillation snapshot takes the same form as Eq. [8] However, the correlation of each multiplied pair 
depends on the extended emission region; the subsequent inclusion of the scintillation ensemble is complicated by 
the different correlations for each scintillation factor: pg ^ pg^i ^ Pg,2- The appropriate correlation and respective 
intensities of the multiplied terms are 
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Also, observe that the mean visibility is weakly diminished by a finite source emission size: 
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3.2. The Relation Between the Emission Region and the Dimensionless Size Parameters js,i 

The precise correspondence between the source dimensions and 7s, ^ requires knowledge of the distribution of scat- 
tering material. The relationship between the correlations {pg, Pg,i} and the observing baseli ne likew ise depends on 
the scattering assumptions. For example, by assuming a square-law phase structure function, iGwinnI (2001, ) derived 
the relations 



7s,i = ( -^k9ia, 



1 - ib^k9y 



exp 



'limef 
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(33) 



Here, D is the characteristic observer-scatterer distance, R is the characteristic source-scatterer distance, k is the 
observing wavenumber, 9i is the angular size of the scattering disk along Si, and at is the standard deviation of 
the (integrated) distribution of source intensity along Si. The scintillation correlations also depend on the baseline 
length |b| and its projections bi along Si. Consequently, for baselines much shorter than the diffractive scale, the 
three correlations are nearly equal. Thus, the dimensionless size parameters 7s, i give the squared size of the source in 
orthogonal directions Si, in units of the magnified diffractive scale. 



3.3. Approximate Evaluation of the PDF of Visibility 

We now apply the i.i.d. approximation to estimate the PDF of visibility, including the effects of an extended emission 
region. Namely, we combine the form of the PDF of visibility within each scintillation snapshot, given by Eq. [S] with 
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Figure 5. Visibility projections for a point source and their modifications (e.g. P(Vr;7s,i) — P(K';7s,i = 0)) corresponding to an elongated 
emission region with 734 = 0.1 and 75,2 = 0.0. The source and background intensities are unity, as are all the pulse amplitude factors Aj. 
The modest averaging, N = 10, reveals the W-shaped signature of emission size observed bv lGwirm et al.l l|2012bl '). 



the distribution of each scintillation random variable: 
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(34) 



Here, we have shifted to polar coordinates for each pair of scintillation gains; e.g. r = {r = |2g|,r' = |zg|,6' = 
arg(zgZg*)}. Eq. [TT]then gives the distribution P{r,r',9) of each triplet. Also, p, /, and /' are defined by Eq. [HH and 
a is then as in Eq. [Sj these variables depend on the integration variables, which account for the scintillation. 

While typical numerical techniqu es for high-dime nsional integrals can evaluate Eq. I34[ they are computationally 
expensive and unenlightening. Gwinn et al.l ()2012b[ ) derived an efficient technique, suitable for N > 20, using the 
Gaussian approximation ( !j2.3.2l) . Their method requires only a one-dimensional integral and a two-dimensional grid 
convolution, which accounts for the signal-dependent noise (see §2.4.2|) . 

However, because we are interested in small values of averaging, especially A^ = 1, we present an alternative. Namely, 
for a small source, anisotropic emission size effects (73^1 7^ 7s, 2) only weakly modify those of an equivalent isotropic 
region. For the analogous effects on intensity, for example, the effects of anisotropy are quadratic in 'ys,i, and so JG12| 
worked in terms of an equivalent isotropic size; 7s = 7s. 1 — 73^2 ■ For our purposes, a one-dimensional emission region is 
advantageous to characterize the dominant effects of size because it obviates the integration over P(zg,2, z' 2), thereby 
reducing the integral of Eq. [34]to six dimensions. 

In addition, by using the results of ^A.2[ we can evaluate projections of the visibility distribution without increasing 
the number of required integrations. For example, the distribution of V^ = Ke{V) is 
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Here, and elsewhere, the subscripts r and i denote the real and imaginary part, respectively. 

Figure [5] shows the effects of an extended emission region on the visibility projections P(T4) and Q(K-), for both 
short and long baselines. 

3.4. The Short Baseline Limit 

The special case of a short baseline vastly simplifies the effects of an extended emission. In this case, the quantities 
in Eq. 1311 onlv depend on the single scintillation "gain" Q: 
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This quantity is the convolution of three independent exponential random variables. The PDF of Q is then (|JG12| . see 
Eq. Al) 
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The PDF of visibility is then given by Eq. [T6l with the substitution G ^f Q, and the projections follow likewise. In 
fact, because p € R, the projections become especially tractable: 
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Figure [5] illustrates how an extended emission region modifies these projections. Increased averaging tends to pronounce 
the effects on the real projection but not the projected imaginary variance because of the d ecreasing var i ance w ith 
averaging. Thus, -P(K) tends to be a more sensitive indicator of size than Q{Vr), as noted bv lGwinn et al.l (|2012b[ ). 

4. SUMMARY 

We have derived the PDF of visibility for a strongly-scintillating source, with particular attention to spectral resolu- 
tion at or near the Nyquist limit. We have incorporated background- and self-noise, source variability, the possibility 
of spatially-extended source emission, and arbitrary temporal averaging. We have also demonstrated that the visibility 
statistics exhibit several remarkable characteristics. For example, the combination of scintillation and self-noise intro- 
duces a broad "skirt" in the distribution of visibility that dominates asymptotic statistics, regardless of the baseline 
or the signal-to-noise. Finally, we have given simplifications of this PDF in various regimes, such as the zero-baseline 
interferometer, as well as results for various projections of the PDF. 

Our results facilitate scintillation studies of pulsars in statistically-delicate regimes and studies of pulsar emission 
regions using interferometry. In particular, our description of Nyquist-limited statistics can provide a sensitive and ro- 
bust detection of an extended emission region, which does not require any assumptions about the nature or distribution 
of the scattering material and can be applied to estimate the emission sizes of individual pulses. 

We thank the referee for a careful reading and for several comments that improved the clarity of the text. We thank 
the U.S. National Science Foundation for financial support for this work (AST-1008865). 



APPENDIX 

A. MATHEMATICAL RESULTS 

A.l. Product of Correlated Complex Caussian Random Variables 

Let wi and W2 be a pair of co rrelated circular complex Gaussian random variables with standard deviations a^ and 
correlation p = [w iw 2) / ijJ 102) ■ iGwinnI (|2001[ ) derived the PDF for the product (i.e. the "visibility") V = wiwj for 
p S R. We generalize his result in two directions: by allowing p Cz C and by accounting for the averaging of N i.i.d. 
visibilities. The first extension is simply a complex rotation of the PDF for real p: 
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Here, we have introduced the scale parameter a = - — ^ — -(Ji(T2 for convenience. The PDF is written with respect to 
the standard complex metric (iRe[K]dIm[T^]. 

The characteristic function of this visibility PDF is then given by 
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where kr and ki are conjugate variables to Re[K] and Im[F], respectively. 
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We can calculate the PDF of the average of N i.i.d. visibilities by inverting the product of their characteristic 
functions: 
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If the averaged visibilities are not statistically isotropic, then we can still obtain a useful reduction using Feynman 
parameters to s ymmetrize the product of characteristic functions. These parameters {si} are defined and applied as 
(|Srednickil[200l : 
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The result of Eq. lA3l can then be applied to the symmetrized integrand. 

Because the visibility is complex, the convolution of N visibilities requires a 2{N — l)-dimensional integral. However, 
to symmetrize the convolution requires a single Feynman parameter for each visibility. The Feynman parameters also 
have an overall (5-function constraint, so the convolution is reduced to an (TV— l)-dimensional integral. If some pairs of 
visibilities are i.i.d., then the dimensionality of the integral can be further reduced. , Gwinn et al. (2012b) applied this 
reduction, in conjunction with the Gaussian approximation (t ^2.3.2p . to evaluate the effects of an extended emission 
region on the PDF of visibility; they thereby reduced the dimensionality of the necessary numerical integration from 
four dimensions to one. 



A. 2. Visibility Projections 

We now calculate the real and imaginary projections of the visibility PDF (Eq. lA3p . We again utilize the characteristic 
function: the conjugate projected variable is set to zero, and the remaining function is inverted with respect to the 
unprojected variable. If the averaged visibilities are i.i.d., we obtain 
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Here, Pr = Re(/9), pi = Im(p), Vi = Re(V^), and fcAr(x) is the modified spherical Bessel function of the second kind 
([Arfken fc WebeHl2005[ ). The imaginary projection P{Vi;a,p,N) follo ws froni the su bstitutions pr ^ A and K- ^ Vi. 
Setting A^ = 1 and p G M recovers the results given in the appendix of lGwimJ (J2001D . 

We also present the projected imaginary variance: Q{Vr) = J V5^P(Vr, Vi)dVi. For this calculation, we again use the 
characteristic function, but take two derivatives with respect to ki and multiply the result by — 1 before zeroing ki and 
inverting with respect to kr'- 
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Additional projections can be performed similarly. 
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